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We report an analysis of decaying into o^tt^/Zt withtj — > Tr^'^Tr^Tr" using a data sample containing 
nearly 320 million r pairs collected with the BaBAR detector at the PEP-II B-Factory. We find no 
evidence for second-class currents and we set an upper limit of 0.69% at 90% confidence level for the 
fraction of second-class currents in this decay mode. 

PACS numbers: 13.35.Dx, 14.60.Fg 



Hadronic weak currents can be classified as either 
first- or second-class depending on the spin J, parity 
P and G-parity G of the final hadronic system In 
the Standard Model, first-class currents (FCC) in t de- 
cays have J^'^ = 0++,0~~,l"' or 1 ^, and are ex- 
pected to dominate. Second-class currents (SCC) have 
J^'^ = 0^ ,0 ^, 1++ or 1 , and are associated with 
a decay constant proportional to the mass difference be- 
tween up and down quarks. Thus they are expected to 



vanish in the limit of perfect isospin symmetry. SCC 
searches have takeivplace extensively in nuclear [3 de- 
cay experiments ||2|,|3|], with no confirmed observations. 
This letter presents a search for SCC in ^ Lun^i'r de- 
cays with uj 'n'^TT~n^ , based on studying the angular 
distributions of final-state particles. 

The decay ijjtt^Vt[^\ is expected to proceed pre- 
dominantly through FCC mediated by the p resonance. 
This decay may also potentially proceed through SCC 
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with J^^ = or 1++. The latter may be mediated 
by 61 (1235) with h^Vr uj-k'Vt. The de- 

cay 6j~ — > WTT" occurs through S- and D-waves |0], as 
compared to a P-wave for FCC. Different ahgnments of 
Lo spin result in different angular distributions of the fi- 
nal state particles. The expected distributions of cos 9^.^, 
F(cos6'(j7r)/ for all possible spin-parity states of the fi- 
nal state particles are listed in Table IB where O^tc is the 
angle between the normal to the lo decay plane and 
the direction of the remaining tt in the lo rest frame. 
The existing measurements of the angular distribution 
in Lon^Vr are consistent with having only the P- 

wave contribution, and the present upper limit is 5.4% 
for the ratio of SCC to FCC contributions at 90% confi- 
dence level (CL) 0,1]. 

TABLE I: Expected angular distributions, F(cos0tj^), for pos- 
sible spin-parity states in the decay — > um'Vr- L is the 
orbital angular momentum. 
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This analysis is based on data recorded by the BaBar 
detector [9] at the PEP-II asymmetric-energy e+e^ stor- 
age rings operated at the SLAC National Accelerator 
Laboratory. The data sample consists of 347.3 fb^^ 
recorded at a center-of-mass energy of 10.58 GeV. With 
a cross section for t pairs of (Jtt ~ (0.919 ± 0.003) nb 
iiol [Till , this data sample contains nearly 320 million 
pairs of t decays. 

The BaBar detector is described in detail in Ref. j^. 
Charged-particle tracks are measured with a 5-layer 
double-sided silicon vertex tracker (SVT) together with 
a 40-layer drift chamber (DCH) inside a 1.5-T supercon- 
ducting solenoid magnet. An electromagnetic calorime- 
ter (EMC) consisting of 6580 CsI(Tl) crystals is used 
for identification of electrons and photons. Charged 
hadrons are identified by a ring-imaging Cherenkov de- 
tector in combination with energy-loss measurements 
(dE/dx) in the SVT and the DCH. An instrumented 
magnetic-flux return (IFR) provides muon identifica- 
tion. 

Monte Carlo simulation is used to estimate the signal 
efficiencies and background contamination. KK2 f [11] is 
used to generate r pairs with the decays of the r leptons 
modeled by Tauola ||12|[. Continuum qq events are sim- 
ulated using JET SET [13]. Final-state radiative effects 
are generated for all decays using Photos [14.1. The de- 
tector response is simulated with GEANT4 |l5|]. and the 
Monte Carlo events are reconstructed in the same man- 
ner as data. 



Since t pairs are produced back-to-back in the e+e~ 
center-of-mass frame, each event is divided into two 
hemispheres according to the thrust axis llT6ll , calcu- 
lated using all reconstructed charged particles. Can- 
didate events in this analysis are required to have a 
"1-3 topology", where one track is in one hemisphere 
(tag hemisphere) and three tracks are in the other hemi- 
sphere (signal hemisphere). Events with four well- 
reconstructed tracks and zero net charge are selected 
for further analysis. The polar angles of all four tracks 
and the neutral clusters used in tt'^ reconstruction are re- 
quired to be within the calorimeter acceptance range. 
Events are rejected if the invariant mass of any pair of 
oppositely charged tracks, assuming electron mass hy- 
potheses, is less than 90 MeV/c^, as these tracks are 
likely to be from photon conversions in the detector ma- 
terial. 

The charged particle found in the tag hemisphere 
must be either an electron or a muon candidate. Elec- 
trons are identified using the ratio of calorimeter en- 
ergy to track momentum (E/p), the shape of the shower 
in the calorimeter, and dE/dx. Muons are identified 
by signals in the IFR and small energy deposits in the 
calorimeter consistent with expectation for a minimum- 
ionizing particle. Charged particles found in the sig- 
nal hemisphere must be identified as pion candidates. 
The tt" candidates are reconstructed from two separate 
EMC clusters with energies above 100 MeV that are 
not associated with charged tracks; these candidates 
are required to have invariant masses between 100 and 
160 MeV/c^. Events are required to have a single 
in the signal hemisphere. The r candidates are recon- 
structed in the signal hemisphere using the three tracks 
and the 7r° candidate. The invariant mass of the r can- 
didate, r7i(7r^7r^7r+7r°), is required to be less than the 
mass of the r lepton. The Monte Carlo simulation pre- 
dicts that 14% of the events remaining after the event 
selection process are r-pair events that do not contain 
a ^ 7r^7r^7r+7r°j/T- decay, and 1.3% are e+e^ qq 
events. 

Each selected event has two tt+tt^tt'^ combinations. 
The LO signal region is defined for masses m(7r+7r^7r°) 
between 760 MeV/c^ and 800 MeV/c^ with mass re- 
gions of width 60 MeV/c^ on each side of the peak used 
as sideband regions for backgroimd studies, as shown 
in Figure [T] For each lo candidate, the angle 6^-,^ is cal- 
culated. The distribution of cos 0^.^ is used for the SCC 
measurement. 

There are three background sources to be consid- 
ered in this analysis: combinatoric background, qq 
events and non-signal t decays. The combinatoric back- 
ground is expected, and confirmed by the simulation, 
to have a distribution of cos Oui-k which is independent 
of r7i(7r+7r^7r°). This allows the sideband regions to be 
used to subtract this background. The number of combi- 
natoric events lying within the signal region is obtained 
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FIG. 1: uj candidate mass spectra for selected events in 
data and background expected from simulation (shaded his- 
togram). The background histogram does not include the non- 
resonant Tr^TT^TT+Tr'Vr decays. The signal (S) and side- 
band (SB) regions are indicated in the figure. 



by fitting the m(7r+7r~7r°) spectrum with a relativistic 
Breit-Wigner convolved with a resolution function for 
the Lo resonance and a polynomial for the combinatoric 
background. The polynomial is integrated over the sig- 
nal region to find the number of continuum events in 
the signal region. 

After subtracting the combinatoric background, ap- 
proximately 0.3% of the remaining events in the signal 
region are expected to be of qq origin, while 4.6% are ex- 
pected to be from non-signal r decays. The background 
from e+e~ qq events that contain lo tt+tt^tt'^ de- 
cays is studied using events with 77i(7r^7r^7r+7r") well 
above the r mass (> 2.1 GcV/c^). In this region, where 
all events are considered to be of qq origin, a compari- 
son of the numbers of w mesons in simulation and data 
is used to scale the simulated qq background before sub- 
tracting from data. 

The dominant non-signal r background, comprising 
99% of the remaining background, is ^ W7r~7r°i^r/ 
where the additional tt" has not been reconstructed. The 
decay uott^tt^Vt has not been well measured 

and is incorrectly modeled in the Monte Carlo. To cor- 
rect for the differences between data and Monte Carlo, 
events with an additional 7r° candidate in the signal 
hemisphere are selected, using the same cuts discussed 
above. Using these events, the Monte Carlo branching 
fraction of Luir^Ti^Vr is corrected by comparing the 

numbers of fitted uj candidates in data and Monte Carlo. 
The fit function used for this is a relativistic Breit-Wigner 
convolved with a resolution function plus a polynomial 
background. The branching fraction obtained using this 
correction technique is found to be consistent with exist- 
ing measurements [6]. To correct the angular distribu- 
tion of ^ LdTi^ TT^Vr, backgrouuds, consisting of com- 
binatorics, events and uj-k^Vt decays (assuming 



the angular distribution corresponding to the dominant 
FCC contribution), are subtracted from the lott^-k^ data 
sample, and the remaining cos 9^^.,^ distribution is used 

distribution in the Monte 



to correct the t' 
Carlo. 

After subtracting backgrounds and applying cos 6*^^- 
dependent efficiency corrections, a binned fit to the re- 
maining cos O^T,: distribution is carried out using 



F(cose'^^^) = iV X 



]-e+-(l-e){l-cos'e^^) 



, (1) 



where iV is a normalization factor and the parameter e is 
the fraction of ^ lott^ Vr decays that proceed through 
sec. In EqllJ only the L = term is used to describe 
the sec contribution since this function gives the most 
conservative estimate of e (i.e. the largest upper limit). 

To help validate the analysis method, the procedures 
were applied to Monte Carlo samples generated to in- 
clude small fractions (1% and 2%) of the second-class 
current process — *■ fei(1235)^i^r lo-k'Vt, as well 
as to a sample containing no SCC contribution. In each 
case the fits to the simulated detector-level, background- 
subtracted angular distributions returned values consis- 
tent with the input fractions of SCC. 

The largest contributions to systematic uncertainties 
on e are scaling and modeling the Monte Carlo back- 
ground. The correction applied to the branching frac- 
tion of ijju^tPvt has an error associated with it, 
determined by the available statistics. This correction 
factor is adjusted by ±1ct to obtain the uncertainty in e 
while the errors associated with correcting the angular 
distribution are folded into the statistical uncertainty. In 
addition, there are r decays that may be present in the fi- 
nal event sample but which are not included in the sim- 
ulation. The largest of these are expected to be 

LuK^Vr, L01T~TT^TT^Vr, and — > WTT^ TT^ TT+Z^^ de- 

cays which, when combined, can add up to 0.2% of the 
final event sample. Since the effect that these decays 
have on the angular distribution is unknown, the ex- 
treme cases are taken to obtain the uncertainty. These 
cases correspond to the decays having either entirely 
1 — cos^ 9i^T^ or entirely cos^ Oujtt distributions. The scal- 
ing of qq events can also affect the measurement of e, 
and the uncertainty is obtained by adjusting the scaling 
factor by ±1(t. These systematic uncertainties are sum- 
marized in Table im 

Subsets of the generic Monte Carlo dataset are used in 
the background studies and in the determination of the 
efficiencies. Therefore, to estimate the sensitivity of the 
analysis without the effect of statistical correlations in 
the Monte Carlo samples, an ensemble of simulated ex- 
periments is used. In this study, angular distributions 
are generated for the signal and sideband regions to 
simulate the statistics available in the data and various 
Monte Carlo samples used in the analysis, with e = in 
the signal Monte Carlo. After subtracting background 
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TABLE II: Summary of systematic uncertainties on e 



xlO^ 



Source 


Uncertainty (a^) 




±0.07% 


non-simulated r decays 


+0.00 (y 
-0.55 


qq scaling 


±0.01% 


Total 


+0.08 0/ 
-0.55 



samples, the angular distribution is corrected for effi- 
ciency and fitted using Eq. [1] The statistical uncertainty 
on e obtained from the fit is 0.63%, which combined with 
the systematic uncertainties leads to an estimated uncer- 



tainty of (Je 



F0.64 
-0.84 



%. 



The angular distribution of the final state particles in 
data is obtained by subtracting estimated backgroimds 
as described above. The remaining distribution is cor- 
rected for efficiency and fitted using Eq. [T] as shown in 
Figure E The fit has x^/dof = 15.4/18, and the fitted 
value of e in the data is (-0.55 ± 0.58(stat)^o;55(syst))%, 
which is consistent with no SCC contribution to — > 
oj-k^Vt decays. 

The upper limit on e is obtained using a Bayesian ap- 
proach [17] with a prior that is flat for e > and zero for 
e < 0. The probability distribution for the value of the 
SCC contribution is a Gaussian with mean e = —0.55% 
and errors <7e =^0 84 taken from the simulation stud- 
ies; however since negative values of e are non-physical, 
only the positive portion of this probability distribution 
is used in the limit calculation. The limits obtained from 
this method are e < 0.68% at 90% CL and e < 0.84% at 
95% CL. 

In summary, a search for second-class currents in 
the decay lott^Vt has been conducted with the 

BaBar detector. No evidence for second-class currents 
is observed, and a 90% confidence level Bayesian up- 
per limit for the fraction of the second-class current in 
T~ uTT^Vr decays is set at 0.68%. For comparison 
with the previous result from CLEO, this is equivalent to 
a ratio of second-class (non-vector) to first-class (vector) 
currents of 0.69%. This limit is an order of magnitude 
lower than the limit set by the CLEO collaboration [8]. 
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FIG. 2: The cos distribution for the data. The curve is the 
result of the fit described in the text. 
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